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A  MATHEMATICAL  MODEL  OF  Till*  PROBABILITY  Oh  PHRFORATION 
OF  Till:  HUMAN  SKULL  BY  A  BALLISTIC  PROJECTILE 


t.  INTRODUCTION 

In  civilian,  as  well  as  military  environments,  there  are  sometimes  hazards  from  ballistic 
fragments.  Because  of  the  extreme  vulnerability  of  the  brain  to  penetrating  injury,  the  difference 
between  a  relatively  harmless  wound  anil  a  lethal  one  often  hinges  on  whether  the  skull  is  perforated 
when  struck  by  one  of  these  projectiles.  In  the  military  application,  an  important  part  of  determining 
the  vulnerability  of  the  soldier  or,  conversely,  the  effectiveness  of  protective  helmets  lies  in  our  ability 
to  predict  the  circumstances  under  which  the  skull  will  withstand  the  impact  of  a  fragment  or  bullet. 
In  this  report  a  mathematical  model  is  developed  which  provides  an  estimate  of  the  probability  of 
perforation  of  the  inner  table  of  the  skull  and  subsequent  penetration  of  the  brain  by  a  cube  or 
sphere.  This  probability  is  considered  t  function  of  the  striking  velocity  and  other  physical  properties 
of  the  projectile  as  well  as  the  thickness  of  the  bone  under  the  impact  area. 

2.  DERIVATION  OF  THE  MODEL 

The  skull  is  composed  of  inner  and  outer  surfaces  of  hard  compact  bone  and  an  area 
between,  called  cancellous  bone,  which  is  less  dense  because  it  contains  a  higher  percentage  of  living 
cells.  The  impacting  projectile,  when  it  perforates,  punches  out  a  hole  in  the  outer  table  roughly  its 
own  size.  The  area  over  which  the  force  is  exerted  spreads  within  the  cancellous  portion,  causing  a 
larger  and  more  irregular  hole  in  the  inner  table.  This  phenomenon  is  often  called  “cratering”  and  is 
used  by  forensic  pathologists  to  tell  exit  from  entrance  wounds  in  the  skull.  This  is  possible  because 
the  larger,  irregular  hole  is  in  the  outer  table  with  an  exit  wound. 

To  simplify  the  mathematical  model,  this  phenomenon  is  idealized  by  the  following 

assumptions: 


1 .  In  the  region  of  the  impact  the  skull  is  a  flat  homogeneous  plate  of  bone  with 
uniform  mechanical  properties. 

2.  When  the  skull  is  perforated  by  a  projectile,  a  truncated  right  circular  cone  of  bone 
is  driven  out  in  front  of  the  fragment. 

3.  The  top  of  the  truncated  cone  has  an  area  equal  to  the  mean  presented  area  of 
the  projectile:  i.e.,  the  projectile  orientation  is  assumed  to  be  random. 

4.  The  velocity  of  sound  in  bone  is  much  greater  than  the  velocity  of  perforation,  so 
compressional  effects  may  be  neglected. 

5.  The  probability  of  perforation  depends  only  on  the  logarithm  of  the  ratio  of  the 
stress  imposed  on  the  skull  by  the  impacting  projectile  to  the  tensile  strength  of  the  bone. 

6.  The  underlying  probability  distribution  governing  perforation  is  well  approximated  by 
the  logistic  distribution  on  the  model  variable,  the  log  stress  ratio  of  assumption  5  (e.g.,  if  the  true 
distribution  were  Gaussian  this  assumption  would  be  valid).  Some  of  the  factors  contributing  to  the 
variance  of  the  probability  distribution  are  variations  in  the  mechanical  properties  of  bone  from  skull 
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ui  skull,  differences  in  the  proportion  of  cortical  to  cancellous  hone,  and  random  deviations  from 
some  of  the  assumptions  listed  above  (for  instance,  deviation  from  conical  shape  in  the  crater),  each 
of  which  is  unique. 


A  schematic  diagram  of  the  penetrating  projectile  is  shown  in  figure  I .  For  purposes  of 
illustration,  the  projectile  is  idealized  as  a  right  circular  cylinder  impacting  end-on.  The  terms  which 
will  he  used  in  developing  the  model,  most  of  which  appear  in  figure  I,  are  defined  in  the  list  of 
symbols.  The  units  used  are  not  SI  hut  are  those  customarily  used  in  wound  ballistics  and  arc 
retained  here  as  a  matter  of  convenience. 


I 

f 


Figure  1.  Idealized  Schematic  Diagram  of 
Cylinder  Perforating  Skull 

The  surface  area  of  the  side  of  the  truncated  cone  is  given  by: 

S  =  ir  (r  +  R)s 

But  s  =  t  sec'P  and  R  =  r  +  t  tan'P,  so 

3  =  ir  (2r  +  t  tan'Pjt  sec^ 


(1) 


C) 
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The  minimum  average  force  exerted  in  pushing  out  the  conical  section  of  skull  may  he  approximated 
by  dividing  the  kinetic  energy  of  the  projectile  by  the  hypothetical  elastic  limit  of  the  bone,  i.c.,  the 
distance  that  the  plug  could  be  moved  before  it  would  break  loose  from  the  surrounding  material. 
Thus, 

force  %  '/imv“/d 


tensile  stress  *=  foree/S 


stress  ratio  =  tensile  stress/tensile  strength  (T) 


_ Vi  mv~/d 

T  tr  (2r  +  t  tan )t  sec'J' 


(3) 


Tiie  quantity  2r  may  be  renlaced  by  the  variable  1)  =  2>JA/ir  to  generalize  to  a  projectile  with  a 
noncircular  presented  area.  When  inserted  into  the  probability  function,  the  multiplicative  constants  of 
equation  3  may  he  absorbed  into  logistic  parameter  a  (see  equation  5  below).  With  these  modifications 
the  stress  ratio  becomes  the  model  variable, 

Vi  mv~  Vi  niv"/t“  * 

*  t~sec^(D/t  +  tan'P)  sec'HD/t  +  tan'l')  <4) 

We  will  use  the  logistic  probability  distribution  function  to  estimate  the  probability  of  penetration  of 
the  skull. 


P  = 


I _ 

I  4.  e-(a  +  blnX) 


(5 1 


where  a  and  b  are  determined  from  data  (see  next  section). 

The  natural  logarithm  (In)  of  X  is  used  instead  of  X  itself  to  equalize  the  variance  over  the  domain  of 
definition.  This  requirement  is  most  easily  explained  by  use  of  an  example:  Suppose  we  are  concerned 
with  the  probability  that  a  certian  projectile  will  perforate  some  material  in  sheets  of  different 
thicknesses  where  the  probability  of  perforation  is  a  function  of  velocity  alone.  On  a  thin  sheet  where 
the  v<5Q,  the  velocity  at  which  the  probability  of  perforation  is  0.5,  is  100  m/s.  we  would  expect  a 
standard  deviation  on  the  order  of  tens  of  meters  per  second.  On  a  much  thicker  sheet,  where  the 
\<;q  is  1000  m/s,  we  no  longer  expect  the  standard  deviation  to  be  measured  in  tens,  but  in  hundreds 
of  meters  per  second.  In  other  words,  we  expect  the  standard  deviation  to  be  roughly  proportional 
to  the  magnitude  of  the  mean.  In  this  case,  the  variance  may  be  approximately  equalized  over  the 
whole  domain  by  dealing  with  the  logarithms  of  the  numbers  rather  than  with  the  numbers  themselves. 


*  The  equation  is  expressed  in  this  manner  as  an  aid  to  plotting  X  at  a  later  stage. 
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The  data  used  were  oblained  from  two  previous  studies  on  skull  penetration.*,**  Since 
the  analysis  on  these  studies  showed  that  the  response  of  dried  skulls  was  significantly  different  from 
that  of  fresh  skullcaps,  only  the  latter  were  used  in  determining  the  parameters  a  and  h  for  the  logistic 
function  of  equation  5.  The  physical  properties  of  the  five  projectiles  used  are  listed  in  the  table. 


I  able.  Projectiles  Studied  (All  Made  of  Steel) 


Projectile 

1 

h— - - - ■ - - - 

Mass 

Mean 

dimension* 

Mean 

presented  area** 

gm 

cm 

cm- 

0.85  grain  sphere 

0.055 

0.238 

0.0445 

2.1  grain  cube 

0.135 

0.265 

0  1050 

4.2  grain  cube 

0.274 

0.333 

0.1660 

id  grain  cube 

1 .029 

0.514 

0.3966 

225  grain  cube 
' 

14.694 

1.236 

2.2933 

*  Diameter  for  sphere:  edge  for  cubes. 
**  Surface  area/4;  see  appendix  A. 


The  ability  of  the  model  to  properly  scale  the  probability  of  perforation  of  skulls  of 
different  thickness  by  projectiles  of  different  mass  was  tested  against  the  v^q  of  the  above  five  pro- 
iectiles.  The  vs0  was  obtained  by  r-eraging  the  five  highest  nonperforating  velocities  and  the  five 
lowest  perforating  velocities  for  each  skull  thickness  for  which  there  were  sufficient  data*  These  data 
are  shown  plotted  in  terms  of  the  model  in  figure  2.  Note  that  the  abseissa  is  the  logarithm  of  the 
numerator  of  the  model  variable  of  equation  4,  while  the  ordinate  is  the  logarithm  of  the  denom¬ 
inator.  A  straight  line  of  slope  1  represents  a  constant  ratio  of  the  numerator  and  denominator  (see 
appendix  B).  When  drawn  through  the  data  in  figure  2,  this  line  allows  a  visual  assessment  of  the 
goodness  of  fit  of  the  data  from  low  energy  to  high  energy.  The  energies  are  mv-$Q  for  all  v^q 
points.  Since  'P  is  a  variable,  it  may  be  adjusted  to  cause  the  sphere  and  cube  data  to  be  on  a 
common  line.  A  T'  of  23°  for  the  sphere  and  30°  for  the  cubes  not  only  superimposed  the  two 
groups  but  tilted  the  individual  data  sets  to  form  the  best  straight  line  of  slope  1 .  Since  the  co¬ 
linearity  of  the  data  with  the  reference  line  is  relatively  insensitive  to  'l/,  the  difference  (7°  between 
spheres  and  cubes)  is  more  significant  than  the  magnitude  of  'P. 


*  Miller.  J.  F  ,  Ashman,  W,  P„  and  Jameson,  J.  W,  Edgcwood  Arsenal  Technical  Report  EATR  4373.  Ballistic  Limits 
of  Skulls  Against  Steel  Cubes.  April  1970. 

**  Miekiewicz.  A.  P.  Chemical  Systems  Laboratory.  Private  communication. 

T  To  obtain  the  maximum  number  of  Vjq  points,  a  few  were  averages  of  4  and  4  instead  of  5  and  5. 


SKULL  BALLISTIC  LIMIT 
X  =  .85  SPHERE  V  =  23  deg 

*  =  2.  I  CUBE  l|l  =  30  deg 

+  =4.2  CUBE 


Alt  .  the  determination  of  all  the  variables  in  equation  4  are  known.  It  is  then 
possible  to  lit  the  dichotomous  data  (perforation  or  nonperforation)  to  the  logistic  function  by  the 
..pproximate  least  squares  method  of  Walker  ami  Duncan  *  In  this  technique  all  data  are  used  to  approx¬ 
imate  the  vjq  lor  each  skull  thickncss/projcctile  combination.  For  each  shot  the  mass,  striking 
velocity.  4'.  skull  thickness,  and  mean  presented  area  are  used  to  calculate  X  from  equation  4.  A 
second  number  is  associated  with  this  value:  I  for  perforation.  0  for  non  perforation.  The  Walker-Duncan 
method  then  iteratively  converges  on  values  of  a  and  h  which  approximately  minimi/e  the  sum  of 
squared  differences  between  the  predicted  probability  of  penetration  (equation  5)  and  the  assigned 
value  (0  to  I)  for  all  shots  by  all  projectiles.  The  resulting  values  of  a  and  b  are: 

a  -  -13.035.  b  =  4.7l>3 

lines  of  slope  I  in  figure  2  represent  constant  values  of  X  and.  therefore,  discrete  levels 
of  probability  of  perforation.  To  find  the  line  of  S09f  probability  we  insert  I*  =  0.5  in  equation  5 

and  solve. 


I 

p  =  - =  o.5 

1  +  e~*a  +  blnX50^ 


or 


e-(a  +  blnXjQ)  =  | 


This  implies 


a  +  blnXjQ  =  0 
1/2  mv~50 

InXcQ  =  In  - - - - — - 

•5U  t“sec4»(D/t  +  tan4') 

where  X^q  is  from  equation  4  with  v  =  vjq. 

Therefore. 


-a/'b 


(b) 


ln[t“sec4'(D/t  +  tan4'))  -  a/b  +  lnl!">  mv^l  0 

It  is  this  line  which  is  drawn  in  figure  2.  The  good  agreement  between  the  line,  derived  from  all  data, 
and  the  plotted  points,  derived  from  a  few  selected  data  points,  lend  some  credence  to  the  validity  of 
using  those  approximate  v$q  points  in  determining  the  value  of  4'.  An  example  of  the  use  of  the 
model  is  given  in  appendix  C. 


*  Walker.  S.  H.,  and  Duncan.  D.  B.  Estimation  of  the  Probability  of  an  Event  as  a  Function  of  Several  Independent 
Variables  Biometrika  J4,  167-179  (1967). 
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4. 


SUMMARY  AND  CONC  LUSIONS 


The  moil  el  ik*  lived  above  is  based  on  the  assumption  that  a  plui:  of  homogeneous  material 
in  the  shape  of  a  truneated  cone  is  driven  out  of  the  skull  by  the  perforating  projectile.  The  area  of 
the  small  end  of  the  truncated  cone  is  assumed  equal  to  the  mean  presented  area  of  the  projectile. 
Under  these  assumptions,  both  spheres  and  cubes  may  be  accommodated  in  the  same  model  by  letting 
tlie  half  angle  of  the  cone  take  on  dilj'cren,  values  for  the  two  shapes.  Experience  with  previous 
models  of  penetration  has  shown  that  clmnkv,  irregular  fragments  usually  behave  in  a  manner  more 
similar  to  cubes  than  spheres.  Therefore,  it  .,  v  .ggested  that  the  cube  half  angle,  'F  =  30°,  be  used 
•dr  fragments  as  weii. 

The  data  used  in  fitting  this  model  were  obtained  from  impacts  on  bare  bone.  In  predicting, 
the  probability  of  penetrating  the  skull  of  live  humans,  this  model  neglects  the  protection  offered  by 
the  scalp  and  (sometimes)  hair.  Thus  the  model  is  conservative  in  the  sense  of  predicting  somewhat 
higher  probabilities  of  penetration  that  would  actually  be  observed  in  live  human  skulls  of  the  same 
thickness. 


Caution  must  be  exercised  when  applying  these  criteria  to  projectiles  larger  than  those  in  the 
data  base  or  of  density  much  different  from  that  of  steel.  There  are  no  data  available  to  test  the 
model  outside  the  range  of  masses  and  presented  areas  of  the  projectiles  of  the  table.  Large  projectiles 
at  lower  velocities  are  particularly  risky  because  of  significant  curvature  of  the  skull  over  large  areas 
and  different  mechanisms  of  skull  fracture  which  occur  under  these  conditions.  It  should  also  be  noted 
that  assumption  3  requires  random  orientation  of  the  impacting  projectile.  This  will  be  valid  for 
projectiles  whose  mean  dimension  is  not  much  larger  than  the  thickness  of  the  skull,  since  in  per- 
lorating  the  skull  the  projectile  will  ultimately  present  its  mean  area,e.g.,by  rotation  during  penetration. 
For  projectiles  whose  mean  dimension  is  greater  than  the  skull  thickness  (for  example,  paper  weights), 
the  projectile  is  unlikely  to  perforate  presenting  its  mean  area,  and  the  effect  of  orientation  of  the 
projectile  will  be  important. 

5.  RECOMMENDATIONS  FOR  FURTHER  WORK 


When  examining  the  skullcaps  perforated  by  cubes,  one  does  n^t  see  flat-sided  holes  with 
neat  square  corners,  but  rounded  holes.  This  indicates  that  the  assumption  of  equality  between  the 
mean  presented  area  of  the  projectile  and  the  area  of  the  hole  in  the  outer  table  should  instead  be 
an  assumption  of  proportionality.  Differences  between  spheres  and  cubes  could  then  be  accounted  for 
by  different  constants  of  proportionality  rather  than  different  half  angles  'F.  Irregular  fragments  would 
have  a  proportionality  constant  somewhere  between  those  of  the  cube  and  sphere.  Unfortunately,  there 
are  no  data  available  for  irregular  fragments  nor  even  sufficient  data  for  the  spheres  and  cubes  to 
allow  an  accurate  determination  of  the  proportionality  constants. 


If  additional  data  are  acquired,  particularly  fragment  data,  it  is  suggested  that  the  following 
model  be  tried: 


Vi  mv^ 

X  =  t(kD  +  t  tan'F) 


where  k  depends  on  the  sh-pe  of  the  projectile  (sphere,  cube,  chunky  fragment,  etc.)  and  'F  is  a 
constant,  independent  of  C!ape. 


APPKNMX  A 


mi:an  i»ri:si:nti:d  arka 

Let  us  first  develop  the  “simple  ease”  of  the  mean  presented  area  of  a  sphere.  We  want  to 
express  this  as  a  fraction  of  its  total  surface  area.  To  obtain  the  surface  area  we  will  integrate  in  polar 
co-ordinates  as  in  figure  A-l.  Notice  that  the  width  of  the  element  of  area  dA  is  r  sin  0  d4>  The  sine 
function  is  necessary  because,  like  the  outer  surface  of  an  orange  segment,  it  must  narrow  to  a  point 
at  the  top  ( 0  =  0).  I  ven  if  the  wedge  is  infinitesimally  thin,  it  still  must  be  shaped  like  a  wedge. 
Another  way  of  looking  at  this  effect  is  to  imagine  0  held  constant;  then  when  we  rotate  around  the 
vertical  (/.)  axis  through  an  angle  ilfk  the  distance  moved  on  the  surface  is  r  sin  0  d$.  The  surface  area 
of  the  sphere  is 


sin<?  d 0  d <f>  =  r^ 


(-cost? )  |*  d$J 


Referring  again  to  figure  A-l,  we  now  wish  to  obtain  the  mean  area  projected  on  the 
plane  PI  at  the  right  of  the  figure,  parallel  to  the  x-z  plane.  The  width  and  height  of  dA  as  projected 
on  this  plane  are  r  sin0  sin^  d <f>  and  r  sin0  d0.  Thus  the  projected  area  is 


Ap 


sin^  sin^f?  d0  d<j> 


(A-2) 


where  the  limits  of  integration  on  <f  are  0  and  ir  because  the  region  from  it  to  2n  is  on  the  far  side 
of  the  sphere  from  the  plane  PI  and  does  not  contribute  to  the  presented  area  required.  Thus. 


The  mean  presented  area  Ap  (as  expected)  is  the  same  as  the  cross-sectional  area  and  is  a  quarter  of 
the  surface  area.  The  angles  6  and  <i>  may  be  interchanged  without  affecting  the  results. 
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Figure  A-l.  Element  of  Area  of  a  Sphere  and 
Its  Projection  on  a  Plane 


We  now  turn  to  the  problem  of  the  mean  presented  area  of  a  polyhedron  with  an  arbitrary 
number  q  of  plane  elements  as  its  surface.  Let  the  area  of  a  typical  element  be  E.  Then  the  projected 
area  of  tile  element  E,  under  all  possible  orientations  of  the  polyhedron  with  respect  to  the  plane  of 
projection  (x-rl.  is  equivalent  to  the  projected  area  of  E  rotated  about  its  center  of  mass  as  in 
figure  A-2.  The  co-ordinate  system  is  centered  at  the  center  of  mass  in  such  a  way  that  the  normal 
n  to  the  area  E  makes  an  angle  9  with  the  z-axis.  The  projection  of  the  normal  onto  the  x-y  plane 
makes  an  angle  <t>  with  the  x-axis.  The  integral  of  the  area  projected  onto  PI  is 


AI  “ 


Jo  Jo 


sin0  sin^  sin0  d0 


(A-4) 


7T 

sin <j>  sin^0  d9  d<6 


7tE 


where  the  first  sin0  is,  as  above,  necessary  to  provide  equal  weight  to  all  possible  orientations.  Because 
1  is  a  finite  area  (not  an  infinitesimal),  the  area  Aj  is  an  integral  projected  area  -  not  a  mean.  To 
obtain  the  mean,  we  must  divide  Aj  by  the  total  solid  angle  subtended  by  all  possible  orientations 
of  the  element  E  (i.e.,  the  solid  angle  swept  out  by  n),  namely  47r.  Thus  the  mean  projected  area 
Ap  =  A,/4t r  =  E/4. 


Appendix  A 
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Figure  A-2.  Element  of  Area  of  a  Polyhedron 
and  Its  Projection  on  a  Plane 


If  there  are  q  elemental  areas  in  the  polyhedron  each  having  a  mean  projected  area  equal 
to  a  quarter  of  its  surface  area,  the  total  mean  projected  area  is  a  quarter  of  the  total  surface  area, 
i.e.. 


q 

=  E 

i=l 


Aj/4 


=  -  V 

4  i=l 


-  4  AS 


f  A-5) 


This  formula  is  valid  for  convex  surfaces  only.  If  there  were  a  concave  portion  on  the 
surface,  it  would  be  shielded  by  another  part  of  the  surface  on  part  of  the  interval  0-ir  of  one  or 
both  angles  of  equation  A-4.  Any  convex  shape  can  be  approximated  to  any  specified  accuracy  by 
making  q  sufficiently  large  and  the  Aj  sufficiently  small.  Thus  the  mean  presented  area  of  a  convex 

1 

solid  is  —  of  its  total  surface  area. 

4 


o 
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APPliNMX  B 


COMMENTS  ON  Till;  SLOP!  OI  A  LOG  -  LOCI  PLOT 

Consider  the  graph  ip  the  figure,  where  the  two  variables  N  and  I)  are  plotted  on 
logarithmic  axes.  Suppose  the  data  indicate  a  straight  line  relationship  on  this  graph. 

Then, 

In  l)  =  p  In  N  +  q 
=  In  N1'  +  q 
In  NP  -  In  D  =  -q 


after  exponentiating. 


NP 

D 


If  p  *  1,  there  is  a  constant  ratio  between  N  and  D. 

If  p  =  2,  there  is  a  constant  ratio  between  and  D. 

If  N  =  !imv-/t~  and  D  =  sec'I'fD/t  +  tan'fr)  the  N/D  =  X  (see  equation  4).  A  line  of  slope  1  in  the 
plot  of  In  N  v  In  D  then  represents  a  constant  ratio  between  Vimv^/t^  and  sec'HD/t  +  tan'I'). 


Figure.  General  Linear  Relationship  Between  In  N  and  In  D 


APPHNDIX  V 


i:xAMPn:  oi  Tin-:  usi:  or  tiii-  matiii-matk  al  modi-l 

Lot  ns  assume  that  we  wish  to  know  the  velocity  corresponding  to  a  0.5  probability  ol 
penetration  lor  a  0.5*gr;un  cube  of  mean  presented  area  0.245  eni“  impacting  a  skull  of  0.5-cm 
thickness.  From  equation  ft  we  obtain 

vqQ  =  |2t“  sec'PfD/t  +  tan*P)exp(-a/b)/in  |  ^ 


where 


I)  =  2  y/sfa  =  2  v'0.245 tn  =  0.5585  cm 


♦  =  30° 


t  =  0.5  cm 

m  =  0.0005  kg 

a  =  -13.035 


b  «  4.793 


Substitution  above  gives  v5q  =  172  m/s.  Suppose  we  wish  to  know  the  probability  of 
penetrating  the  skull  if  the  cube  of  our  example  impacts  at  200  m/s.  From  equation  4  we  obtain 


X  = 


Vi  mv^ _ 

t^sec'PfD/t  +  tamP) 


20.44 


Then 


1 

p  - - - - =  o  807 

1  +  exp(-a  -  blnX) 

Tlie  following  table  can  be  constructed  using  similar  calculations,  assuming  the  same  fragment  and 
skull  thickness: 


Velocity  m/s 


Probability  of  penetration 


150 

175 

200 

225 


0.209 

0.537 

0.807 

0.973 


The  value  of  probability  of  penetration  above  does  not  indicate  the  proportion  of  the  skull  thickness 
penetrated  but  rather  the  proportion  of  hits  which  would  be  expected  to  punch  out  a  plug  of  bone 


from  the  skull  given  a  large  number  of  hits.  Thus,  in  the  example  above,  at  v  =  200  m/s,  80  7%  of 
all  shots  bred  would  be  expected  to  penetrate  for  that  fragment /skull  thickness  combination,  while 
l -  would  not.  No  attempt  is  made  to  .-redict  the  outcome  of  individual  shots 
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0 I STRIBUT I  ON  LIST  II 


Names  Copies 


CHEMICAL  SYSTEMS  LABORATORY 

ATTN: 

DRDAR-CLF 

1 

ATTN: 

DRDAR-CLC-B 

1 

ATTN: 

DRDAP-CLC-C 

1 

ATTN: 

DRDAR-CLC-E 

1 

ATTN: 

DRDAR-C'.J-R 

3 

ATTN: 

DRDAR-CLJ-L 

3 

ATTN  • 

DRDAR-J-M 

1 

ATTN: 

DRDAR-CL J-P 

1 

ATTN: 

DRDAR-CLT 

1 

ATTN: 

DRDAR-CLT-E 

1 

ATTN: 

DRDAR-CLT-P 

1 

ATTN: 

DkDAR-CLN 

1 

ATTN: 

0RDAR-CLW 

1 

A  T  T  N  : 

DRDAR-CLW-C 

! 

ATTN: 

DRDAR-CLW-P 

1 

ATTN: 

DRDAR-CLW-E 

1 

ATTN: 

DRDAR-CLB-C 

1 

ATTN: 

0RDAR-CLB-P 

1 

ATTN: 

DRDAR-CLB-PA 

1 

ATTN: 

DRDAR-CLB-PO 

1 

ATTN: 

DRDAR-CLB-R 

1 

ATTN: 

DRDAR-CLB-T 

1 

ATTN: 

DRDAR-CLB-TE 

1 

ATTN: 

DRDAR-CL Y-A 

1 

ATTN: 

DRDAR-CLY-R 

6 

ATTN: 

0RDAR-CLR-  1 

1 

COPIES 

FOR  AUTHOR(S): 

Research  Division 

12 

DEPARTMENT  OF  DEFENSE 

Defense  Technical  Information  Center 
ATTN:  OT I C-DDA-2  12 

Cameron  Sta+ion,  Building  5 
Alexandria,  VA  22314 

D 1  rector 

Defense  Intelligence  Agency 

ATTN:  DB-4GI  I 

Washington,  DC  20301 

Special  Agent  in  Charge 

ARO,  902d  Military  Intelligence  GP 

ATTN:  I AGPA-A-AN  1 

Aberdeen  Proving  Ground,  MO  21003 


Names  Copies 


Commander 
SED,  HQ,  INSCOM 

ATTN:  I  RFM-SF.D  (Mr.  Joubert)  1 

Fort  Meade,  MD  20755 

DEPARTMENT  OF  THE  ARMY 

HQDA  (DAM0-NCC)  I 

WASH  DC  20310 

Federal  Emergency  Management  Agency 
Office  of  Mitigation  and  Research 

ATTN:  David  W,  8ensen  I 

Washington,  DC  20472 

Deputy  Chief  of  Staff  t  Research, 
Development  &  Acquisition 
ATTN:  DAMA-CSS-C  1 

ATTN:  DAMA-ARZ-D  1 

Washington,  DC  20310 

Department  of  the  Army 
Headquarters,  Sixth  US  Army 
ATTN:  AFKC-OP-NBC  I 

Presidio  of  San  Francisco,  CA  94129 

US  Army  Research  and  Standardization 
Group  (Europe) 

ATTN:  DRXSN-E-SC  1 

Box  65,  FP0  New  York  09510 

HQDA  (DAMI-FIT)  1 

WASH,  DC  20310 

Commander 

HQ  7th  Medical  Command 

ATTN:  AEFPM  1 

AP0  New  York  09403 

Commander 
DARC0M,  STITEUR 

ATTN:  DRXST-ST1  1 

Box  48,  A0  New  York  09710 

Commander 

USAOTEA 

ATTN:  CSTE-ZX  1 

5600  Columbia  Pike 
Falls  Church,  VA  22041 
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Commander 

US  Army  Science  &  Technology  Center- 
Far  East  Office 

ATTN:  MAJ  Borges  1 

APO  San  Francisco  96328 


Commander 

2d  Infantry  Division 

ATTN:  EAIDCOM  1 

APO  San  Francisco  96224 

Commander 

5tn  Infantry  Division  (Mech) 

ATTN:  Division  Chemical  Officer  1 

Fort  Polk,  LA  71459 

OFFICE  OF  THE  SURGEON  GENERAL 

Comma  nder 

US  Army  Medical  Bioengineering  Research 
and  Development  Laboratory 
ATTN:  SGRD-UBD-AL  1 

Fort  Detrlck,  Bldg  568 
Frederick,  MD  21701 

Headquarters 

US  Army  Medical  Research  and 
Development  Command 

ATTN:  SGRD-PL  1 

Fort  Detrlck,  MD  21701 

Commander 

USA  Biomedical  Laboratory 

ATTN:  SGRD-UV-L  1 

Aberdeen  Proving  Ground,  MD  21010 

US  ARMY  HEALTH  SERVICE  COMMAND 


US  ARMY  MATERIEL  DEVELOPMENT  AND 
READINESS  COMMAND 

Commander 

US  Army  Materiel  Development  and 
Readiness  Command 

ATTN:  DRCLDC  1 

ATTN:  DRCSF-P  1 

5001  Elsenhower  Ave 
Alexandria,  VA  22333 

Project  Manager  Smoke/Obscurants 

ATTN:  DRCPM-SMK  1 

Aberdeen  Proving  Ground,  MD  21005 

Director 

Human  Engineering  Laboratory 

ATTN:  DRXHE-SP  (CB  Defense  Team)  1 

Aberdeen  Proving  Ground,  MD  21005 

Commander 

US  Army  Foreign  Science  &  Technology 
Center 

ATTN:  DRXST-MT3  ' 

220  Seventh  St.,  NE 
Charlottesville,  VA  22901 

Dl  rector 

US  Army  Materiel  Systems  Analysis  Activity 
ATTN:  DRXSY-MP  1 

ATTN:  DRXSY-T  (Mr.  Metz)  2 

Aberdeen  Proving  Ground,  MD  21005 

Commander 

US  Army  Missile  Command 
Redstone  Scientific  Information  Center 
ATTN:  DRSMI-RPR  (Documents)  1 

Redstone  Arsenal,  AL  35809 


Super  I ntendent 
Academy  of  Health  Sciences 
US  Army 

ATTN:  HSA-CDH 

ATTN:  HSA-IPM 

Fort  Sam  Houston,  TX  78234 

US  Army  Natick  Research  and 
Development  Command 


ATTN: 

DRDNA-VC 

ATTN: 

DRDNA-VCC 

ATTN: 

DRDNA-VM 

ATTN: 

drdna-vr 

ATTN: 

0RDNA-VT 

Nat  1 ck 

,  MA  01760 

Dl rector 

DARC0M  Field  Safety  Activity 
ATTN:  DRXOS-C  1 

1  Charlestown,  IN  47 1 1 1 
1 

Commander 
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US  ARMY  ARMAMENT  RESEARCH  AND 
DEVELOPMENT  COMMAND 

Comma  ndor 

US  Army  Armament  Research  and 
Development  Command 
A 1 ( N :  DROAR-LCA-L 
ATTN:  DRDAR-LCE 
A!  IN:  DRDAR-I.Ct -C 
ATTN:  ORDAR-I.CM-E 

ATTN:  DRDAR-LCM-SA 
ATTN;  I1R0AR-LCU 
AHN:  DROAR-LCU-CE 
MFN:  ORDAR-PMA  (G.R,  Sacco) 

ATTN:  ORDAR-SCM 
AHN:  Dnl.'.AR-SCP 
ATTN:  HRD AR- 5CA-W 

ATTN:  PRDAR-SER 
ATTN.  DRDAR-TSS 
A'  IN;  DRDAR-CAWS-AM 
A f  r N :  DROAR-CAWS-S I 
Do  vo  ' ,  NJ  07601 

D i rector 

Ballistic  Research  Laboratory 

ARRAOCOM 

ATTN:  DRDAR-TSB-S 

Aberdeen  Proving  Ground,  MD  21005 

D I  rector 

Benet  Weapons  Laboratory 
ATTN:  DRDAR-i  CB-TL 

Water  / 1  I  et,  l,Y  12  189 

Commander 

USA  Technical  Detachment 
US  Naval  EOD  Facility 
Indian  Head,  MD  20640 

US  ARMY  ARMAMENT  MATERIAL  READINESS 
COMMAND 

Commander 

US  Army  Ar nament  Merer  lei 
Readiness  Command 
ATTN:  DRSAR-ASN 
ATTN:  DRSAR-IRC 
ATTN:  DRSAR-ISE 
ATTN:  DRSAR-LEP-L 
ATTN :  DRSAR-PDM 
ATTN:  DRSAR-PE 
ATTN:  DRSAR-SF 
ATTN:  DRSAR-SR 
Rock  Island,  IL  61299 


Commander 
USA  ARRCOM 
ATTN:  DRSAR-MAS-C 
ATTN:  DRSAR-MAD-E 
ATTN:  SARTE 

1  Aberdeen  Proving  Ground,  MD  21010 

1 

)  Commander 

I  US  Army  Dugway  Proving  Ground 

1  ATTN:  Technical  Library  (Docu  Sect) 

I  Dugway,  U1  84022 

1 

l  US  ARMY  TRAINING  4  DOCTRINE  COMMAND 

1 

t  Commandant 

1  US  Army  Infantry  School 

1  ATTN:  NBC  Division 

2  Fort  Banning,  GA  31905 
t 

1  Commandant 

US  Army  Missile  &  Munitions  Center 
and  School 
ATTN:  ATSK-CD-MD 
ATTN:  ATSK-DT-MU-E0D 
Restone  Arsenal,  AL  35809 

1 

Commander 

US  Army  Logistics  Center 
ATTN:  ATCL-MG 
Fort  Lee,  VA  23801 

1 

Commandant 
USAMP&CS/TC&FM 
ATTN:  ATZN-CM-CDM 
1  Fort  McClellan,  AL  36205 

Commander 

US  Arny  Infantry  Center 
ATTN:  ATSH-CD-fS-C 
Fort  Banning,  GA  1905 

Commander 

US  Army  Infantry  Center 
ATTN:  ATZB-DPT-P0-NBC 
’  Fort  Banning,  GA  31905 

1  Commando. 

1  USA  Training  and  Doctrine  Command 

1  ATTN:  ATCD-Z 

1  Fort  Monroe,  VA  23651 


Commander 

USA  Combined  Arms  Center  and 
Fort  Leavenworth 

ATTN:  ATZL-CA-COG  I 

ATTN:  ATZL-CAM- I M  1 

Fort  Leavenworth,  K$  66027 

Commander 

US  Army  TRADOC  System  Analysis  Activity 
ATTN:  ATAA-SL  1 

Whito  Sands  Missile  Range,  NM  86002 

US  ARMY  TEST  &  EVALUATION  COMMAND 

Comma  nder 

US  Army  Test  A  Evaluation  Command 
ATTN:  DKSTE-CM-F  I 

ATTN:  DRSTE-CT-T  I 

Aberdeen  Proving  Ground,  MO  21005 

DEPARTMENT  of  the  navy 

Chief  of  Naval  Research 

ATTN:  Code  443  I 

800  N.  Quincy  Street 
Arlington,  VA  22217 

Commander 

Naval  Explosive  Ordnance  Disposal  Facility 
ATTN:  Army  Chemical  Officer  (Code  AC-3)  1 

Indian  Head,  MD  20640 

Commander 

Naval  Surface  Weapons  Center 
Wh i te  Oak  S I te 

ATTN:  Library  1-321  1 

SI  I ver  Spring,  MD  20910 

Chief,  Bureau  of  Medicine  &  Surgery 
Department  of  the  Navy 

ATTN:  MED  3C33  1 

Washington,  DC  20372 

Commander 

Naval  Weapons  Center 

ATTN:  Technical  Library  (Code  343)  1 

China  Lake,  CA  93555 

Commander  Officer 

Naval  Weapons  Suport  Center 

ATTN:  Code  5042  (Dr.  8.E.  Doude)  1 

Crane,  IN  47522 


US  MARINE  CORPS 

Commandant 

HQ,  US  Marine  Corps 

ATTN:  Code  LMW  I 

Washington,  DC  20380 

Director,  Development  Center 
Marine  Corps  Development  and 
Education  Command 

ATTN:  Fire  Power  Division  1 

Quant  I co,  VA  22134 

DEPARTMENT  OF  THE  AIK  FORCE 

HQ  Foreign  Technology  Division  (AFSC) 

ATTN:  TQTR  I 

Wrlght-Patterson  AFB,  OH  45433 


Commander 

Aeronautical  Systems  Division 

ATTN:  ASD/AELD  1 

ATTN:  ASD/AESD 

Wrlght-Patterson  AFB,  OH  45433 

HQ  AFLC/LOWMM  1 

Wrlght-Patterson  AFB,  OH  45433 

HQ,  AFSC/SDNE  1 

Andrews  AFB,  MD  20334 

HQ  AMD/RD 

ATTN:  Chemical  Defense  OPR  1 

Brooks  AFB,  TX  78235 

HQ  AFISC/SEV  1 

Norton  AFB,  CA  92409 

N0RAD  Combat  Operations  Center 

ATTN:  D0UN  1 

Cheyenne  Mtn  Complex,  CO  80914 


Air  Force  Aerospace  Medical  Research 
Laboratory 

ATTN:  AFAMRL/HE  (  Dr.  C.R.  Replogle)  1 

Wrlght-Patterson  AFB,  OH  45433 

USAF  SAM/RZW  1 

Brooks  AFB,  TX  78235 

HQ  AFTEC/SGB  1 

Klrt land  AFB,  NM  87117 
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OUTSIDE  AGENCIES 


ADDITIONAL  ADDRESSEES 


Dattelle,  Columbus  Laboratories 
ATTN:  TACTEC 
5  0  5  King  A  .  .  ;hj  o 
Columbus,  "H  43201 

Toxicology  Information  Confer,  WG  1008 
Niii-tnal  Research  Council 
2  1 0  1  Constitution  Ave,  ,  NW 
Washington,  DC  20418 

rS  Public  Health  Service 
Canter  for  Disease  Control 
A  Tt  n  .  lewis  Webb,  Jr. 

Du  1  I  u  i  II.;  4  ,  Room  232 
■'!  I  duta,  G.-  30333 

.  1  •  'C  ♦  'r 

Coutoi  u*  lllgence  Agency 

ATTN;  0KC/DD/S&T 
Wasnington,  DC  20503 


Common dor 

1  IIS  Army  1  n  v  I  ronmenta  I  Hygiene  Agency 

ATTN:  librarian.  Bldg  2100 

Aberdeen  Proving  Ground,  MD  21010 

St  I mson  Library  (Documents) 

1  Academy  of  Health  Sciences 

Bldg.  2040 

Tort  Sam  Houston,  TX  78234 
PI  rector 

Ballistic  Research  Laboratory 
1  ARRADCOM 

ATTN:  DROAR-BL  V 

Aberdoon  Proving  Ground,  MD  21005 


1 


